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I.  INTRODUCTION 


Acoustic  scattering  from  the  ocean  floor  has  been  an  item  of  considerable 
interest  in  recent  years.  Since  there  is  no  general  agreement  on  what  the  basic 
scattering  mechanisms  are,  a  lot  of  models  have  been  developed,  based  on 
several  different  scattering  mechanisms.  Some  models'*  >2  assume  that  the 
scattering  occurs  at  the  fluid/sediment  interface.  Other  models^-^  consider 
scattering  from  within  the  sediment  volume  below  the  interface.  Volume 
scatterers  that  have  been  considered  include  variations  in  the  refractive  index  of 
the  sediment  caused  by  local  variations  of  the  sound  speed  or  density,  the 
sediment  grains,  fluid  pockets  between  grains  and,  most  recently,  trapped  gas 
bubbles.  It  is  reasonable  that,  when  gas  is  present  in  sufficient  quantity,  the  gas 
bubble  scattering  might  dominate,  since  gas  bubbles  in  resonance  are  very 
strong  scatterers  of  sound,  with  scattering  cross  sections  that  are  typically 
1000  times  their  geometric  cross  section.®  The  model  of  Boyle  and  Chotiros® 
suggests  that  very  small  amounts  of  gas  (fluid  gas  fractions  of  1 0'®  or  less)  are 
sufficient  to  dominate  the  other  scattering  mechanisms. 

A  resulting  difficulty  in  modeling  acoustic  backscatter  is  that  the 
backscattering  strength  of  a  gassy  sediment  depends  sensitively  on  the  fluid  gas 
fraction.  There  is  at  present  no  practical  method  of  measuring,  in  situ,  small 
amounts  of  gas. 

Actual  sediments  may  contain  more  than  one  type  of  scatterer,  with 
acoustic  returns  that  appear  quite  similar.  For  example,  a  hard  rock  surface  may 
look  very  much  like  a  pocket  of  trapped  gas.  This  problem  might  seriously  affect 
any  attempt  at  remote  acoustic  classification  of  sediment.  Therefore  finding  a 
way  of  distinguishing  gas  bubble  acoustic  returns  from  other  types  is  critical. 

One  way  in  which  gas  bubbles  differ  from  other  kinds  of  scatterers  is  in 
their  nonlinear  scattering  properties.  Water  has  a  nonlinearity  parameter  B/A  of 
about  5-6  compared  with  a  value  of  8-12  for  water  saturated  sand.'*®  The  value 
for  the  gas  inside  a  bubble  is  much  higher.*  1 

One  possible  way  of  identifying  bubble  returns  would  be  to  insonify  the 
sediment  with  a  parametric  signal  and  measure  the  scattered  sound  at  the 
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difference  frequency.  This  procedure  has  been  successfully  employed  to  detect 
bubbles  in  the  water  colurnn.^^ 

The  problem  of  identifying  gas  bubbles  trapped  within  sediment  pores  is 
considerably  more  complicated  than  identifying  bubbles  in  the  water  column. 
Recent  experiments'' ^.1 4  suggest  that  sandy  sediments,  which  are  common  in 
coastal  areas,  are  best  modeled  acoustically  via  Biot's^ poroelastic  theory. 
This  type  of  treatment  differs  from  others  in  that  it  models  a  saturated  sediment 
as  a  two-phase  medium,  consisting  of  a  semi-rigid  structure  of  sediment  grains 
through  which  the  pore  fluid  is  allowed  to  flow.  An  outstanding  feature  of  Biot's 
theory  is  that  two  compressional  waves  are  predicted,  in  addition  to  the  shear 
wave.  In  the  so-called  Biot  fast  wave,  the  pore  fluid  and  grain  structure  move 
approximately  in  phase.  This  wave  is  analogous  to  the  compressional  wave  that 
would  exist  in  an  equivalent  elastic  medium.  The  other  type  of  compressional 
wave,  called  the  Biot  slow  wave,  consists  of  the  solid  and  fluid  parts  of  the 
medium  moving  out  of  phase.  The  experiments  suggest  that  both  types  of  waves 
can  be  significant.  For  trapped  bubble  scattering  calculations,  shear  waves  can 
be  neglected  since  they  couple  only  weakly  into  the  pore  fluid  that  contains  the 
scattering  bubbles. 

As  a  bubble  expands,  it  pushes  against  the  surrounding  pore  fluid.  Partial 
confinement  within  sediment  pores  will  cause  the  fluid  to  behave  as  if  it  had  an 
effective  density  different  from  its  actual  density.  Since  the  mechanics  of  fluid 
confinement  are  different  for  Biot  fast  and  slow  waves,  the  pore  fluid  can  be 
expected  to  have  two  different  effective  densities,  depending  on  whether  fast  or 
slow  waves  are  propagating.  Similarly,  the  fluid  acoustic  impedance  and  the 
bubble's  resonance  frequency  and  damping  constant  are  specific  for  each  of  the 
Biot  waves.  Furthermore,  whereas  for  a  fluid  model  these  quantities  can 
generally  be  treated  as  constants,  in  the  Biot  case  they  are  functions  of 
frequency. 

In  this  report,  an  expression  for  the  difference  frequency  scattering  cross 
section  from  a  bubbly  sediment  is  developed.  It  separately  models  scattering 
from  Biot  fast  waves  and  Biot  slow  waves  by  modeling  the  medium  as  a 
superposition  of  two  acoustic  fluids.  One  fluid  supports  the  fast  wave  while  the 
other  supports  the  slow  wave.  This  is  attained  by  assigning  an  appropriate 
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effective  density  and  acoustic  impedance  to  each  fluid.  The  densities  and 
impedances  are  obtained  from  the  Biot  model  of  Stern,  Bedford,  and  Millwater.''^ 
Shear  waves  are  neglected  since  they  exist  primarily  in  the  solid  part  of  the 
medium  and  do  not  couple  strongly  into  the  pore  fluid  that  contains  the  scattering 
bubbles. 

In  Section  2,  an  expression  for  the  difference  frequency  scattering  cross 
section  from  a  single  bubble  in  an  acoustic  fluid  is  derived.  The  bubble's 
resonance  frequency  and  damping  constant,  as  well  as  the  surrounding  fluid's 
density  and  acoustic  impedance,  are  allowed  to  vary  with  frequency.  In  Section  3 
the  scattering  cross  section  from  a  distribution  of  bubble  sizes  is  computed. 
Section  4  contains  a  calculation  of  the  effective  backscattering  strength  of  the 
interface  above  a  bubbly  Biot  medium.  Conclusions  are  discussed  in  Section  5. 
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2.  CALCULATION  OF  THE  DIFFERENCE  FREQUENCY  SCATTERING 
CROSS  SECTION  FROM  A  SINGLE  BUBBLE 


Ostrovsky  and  Sutin^®  developed  an  expression  for  the  difference 
frequency  scattering  cross  section  from  bubbles  in  an  acoustic  fluid.  The 
following  development  follows  their  technique,  allowing  for  the  frequency 
dependence  of  the  pore  fluid's  effective  density  and  acoustic  impedance  and  the 
bubble's  damping  constant  and  resonance  frequency. 

In  Appendix  A,  Zabolotskaya  and  Soluyan's^®  equation  of  motion  for  a 
bubble's  surface  is  derived: 


V  +  a)o2V-aV2-p(2W+\/2)  +  co5\/  =  -eP  ,  (2.1) 

where  V  is  a  perturbation  in  the  bubble's  volume,  P  is  the  acoustic  pressure 
incident  upon  the  bubble,  co  is  the  frequency  of  volume  oscillations,  and  coo,  a,  p, 
8,  and  e  are  real  quantities  that  are  allowed  to  vary  with  frequency.  The  bubble's 
radius  is  assumed  small  in  comparison  with  the  acoustic  wavelength.  A  solution 
is  needed,  given  a  bi-harmonic  acoustic  pressure  incident  upon  the  bubble; 

P  =  PiCOS(C0it  +  (P^)  +  P2COS(tD2t  +  (P2)  ,  (2.2) 

where  Pi  and  P2  are  amplitudes  and  (pi  and  (p2  are  corresponding  phases  of  two 
superimposed  incident  pressure  signals  with  frequencies  cdi  and  a)2-  As  an 
ansatz,  let  V  be  of  the  form 

V  =  Vicos(coit  +  ^1)  +  V2cos(o)2t  +  +  VnCOs(Qt  +  ^n)  +  ^  ,  (2.3) 

where  Q  =coi-co2  is  the  difference  frequency.  Vj  and  are  amplitudes  and  phases 
of  superimposed  components  of  V.  "'F"  includes  all  other  terms,  including  those 
involving  the  sum  frequency  and  higher  harmonics  of  the  driving  frequencies  coi 
and  (02.  To  solve  for  Vi  and  V2,  it  is  only  necessary  to  consider  the  linear  part 
ofEq.  (2.1): 


5 


V  +  (Oq^  +  o)8V  =  -eP 


(2.4) 


where  coo,  8,  and  e  are  functions  of  the  frequency  to  of  V.  The  linear  part  of  V  is 
given  by 


V|in  =  ViCOs(a)it  +  +  V2COs(a)2t  +  y  .  (2.5) 


Upon  substituting  Eqs.  (2.5)  and  (2.2)  into  Eq.  (2.4), 


(cooi^  -  o)i2j\/^cos(Q)it  +  ^i)  +  5iCDi^iSin(ci)it  +  ^i) 

+{o)02^  -  Q)2^)V2COS(o)2t  +  y  +  820)2^221^(0)2^  +  yj 


-eiPiCos(o)it  +  cpi) 
-e2P2COS(o)2t  +  (P2) 


(2.6) 


If  Eq.  (2.6)  is  valid  for  all  possible  values  of  time  t,  the  o)i  and  0)2  components 
must  satisfy  it  separately.  Hence, 


(o)oi^  -  o)i2)ViCOs(o)it  +  ^1)  +  8iO)i2ViSin(Q)it  +  ^1)  =  -eiPiCOs(o)it  +  (p^)  (2.7) 

(0)02^-O)2^)V2COS(0)2t  +  y +  820)2^2Sin(0)2t  +  y  =  -£2P2C0S(0)2t  +  (p2)  .  (2.8) 


By  equating  magnitudes  on  either  side  of  Eqs.  (2.7)  and  (2.8),  the  amplitudes  Vi 
and  V2  are  obtained: 


Vi-  - !l5 _ 

+  51^0)/ 

(2.9) 

V2=  _ 

V  (®02^  “  ^2^  +  ^2^2^ 

(2.10) 

In  order  to  find  the  difference  frequency  component  Vq  of  the  bubble 
volume  perturbation  V,  it  is  necessary  to  consider  the  complete  nonlinear 
equation,  Eq.  (2.1).  Each  term  in  this  equation  will  have  a  contribution  to  Vq. 
The  first  nonlinear  term  in  Eq.  (2.1)  is  (-aV2).  To  find  the  (o)i-a)2)  component  of 
V2  ,  apply  expression  (2.3): 
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(2.11) 


=  Vi  ^COS^fcDit  +  ^i)  +  2V1V2  cos  (cDit  +  ^1)  cos  (o)2t  +  ^2) 

+  V2^cos^(o)2t  +  ^2)  +  higher  order  terms 

By  invoking  the  trigonometric  identity 

cos(a)cos(b)  =  ^|cos(a+b)  +  cos(a-b) j ,  (2. 1 2) 

the  difference  frequency  component  of  can  be  computed: 

(v%  =  V,V2COs[nt  + (4,-^2)].  (2.13) 

In  similar  fashion,  the  difference  frequency  components  of  the  other  nonlinear 
terms  in  Eq.  (2.1)  can  be  computed: 


2  ^V,V^C08[QH-(^,-^2)j  (2.14) 

(V^)Q  =  m,ci)2V,V2COs[nt  +  (4,-52)j  .  (2.15) 

The  linear  terms  of  Eq.  (2.1)  also  have  Q  components: 

(V)Q  =  VaCOs(Qt  +  ^3)  (2.16) 

(V)(,  =  £2VnSin(ni  +  43)  (2.17) 

('5')Q  =  -a%cos(£2t  +  |3)  .  (2.18) 


The  total  Q  component  of  Eq.  (2.1)  is  the  sum  of  contributions  from  each  of  its 
terms.  Therefore,  by  substituting  Eqs.  (2.13  -  2.18)  into  Eq.  (2.1)  and  noting  that 
the  incident  pressure  P  has  no  difference  frequency  component. 


{-  Q^Vj2Cos(iit  +  ^3) 

+  COo£2%COs(C2t  +  ^3) 

-  ttQ  Vi  V2  cosj^iit  +  (^1-^2)] 

-  Pn  (-t0i-®2)ViV2  cos[i2t  +  (^1-^2)] 

~  pQ  CO102V1V2  COsj^Qt  +  {^1-^2)] 

+  6Qa%sin(Qt  +  ^3) }  =  0  .  (2.19) 

where  cooq,  a^,  |3q,  and  5n  are  the  values  of  coq,  a.  P,  and  8  at  the  difference 
frequency  Rearranging, 

(®OQ^^^)yQCOs(i2t  +  ^3)  +  5Qi2^VQSin(iit  +  ^3) 

=  («£2  -  Pn(o)f+a)|-coiCo^)ViV2  co^ilt  +  •  (2.20) 

Expression  (2.20)  is  an  equality  of  harmonic  signals.  For  two  harmonic  signals  to 
be  equal,  their  magnitudes  and  phases  must  be  equal.  By  equating  the 
magnitudes  on  either  side  of  Eq.  (2.20), 

I  \l Vq I  =  I  {aj2  -  p Ja)f+o)|-a)iCt)2))Vi V2 1  ;  (2.21) 

hence  the  amplitude  of  the  difference  frequency  oscillations  of  the  bubble  volume 
V  is 


|Vd 


(«n  -  PnK+coi-coia)2))ViV2 


(2.22) 
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Upon  inserting  expressions  (2.9)  and  (2.10)  for  Vi  and  V2, 


r 

ttn-Pn 

COl+COl-COiCOzj 

l^i^aPiP 2 

*ny{o)„'-o),f+8, 

(2.23) 

This  expression  reduces  to  that  of  Zabolotskaya  and  Sutin^o  for  the  special  case 
where  coo,  a,  p,  and  6  are  independent  of  frequency. 

As  the  bubble's  volume  oscillates  at  frequency  Q,  a  pressure  PQ(r)  is 
generated  in  the  surrounding  fluid  at  distance  r  from  the  bubble's  center.  In 
Appendix  B,  an  equation  is  derived  that  relates  this  pressure  to  the  bubble's 
volume  perturbation  Vq: 


|Po(r)l  = 


47tr 


Upon  substituting  Eq.  (2.23)  for  IVqI  in  Eq.  (2.24), 


(2.24) 


|Po(r)l  = 


Q  -  Pn(cof+coi-coiC02))eie2Pi  P2 
AnrsJ  sj  {coqi^  -  0)i2)^+5i^a)/  sj  (coo2^  - 


(2.25) 

The  scattering  cross  section  of  the  bubble  is  the  time  average  of  the  scattered 
power  divided  by  the  incident  intensity: 


o  =  47tr^ 


Pi  V1+P2V2 


(2.26) 
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where  Pi  and  P2  are  the  acoustic  pressure  amplitudes  carried  by  the  incident 
waves  upon  the  bubble,  defined  in  Eq.  (2.2).  The  quantities  vi  and  V2  are  the 
corresponding  fluid  velocities.  Pn(r)  and  V£i(r)  are  the  pressure  and  velocity 
amplitudes  carried  by  the  scattered  wave  from  the  bubble  at  distance  r. 


In  terms  of  local  acoustic  impedances  Zj=dPi/dVi, 


a  =  47ur^ 


(Pn(r))^ 

L-\  L2 


(2.27) 


Substitution  of  Eq.  (2.25)  for  Pn(r)  in  Eq.  (2.27)  yields 


_ QV^«n-Pn(ci)f+co|-coiC02))^efeiPfP| _ 

47i[(tOon^Q^)%5QQ'^J(o)oi^cof)^+5fo)|J(a)o2^^i)^+52C02|pf^+P2^ 


.(2.28) 


Upon  substituting  the  following  expressions,  defined  in  Appendix  A, 


47lReq 

=  3Po(Y+ 1)^00  > 


(2.29) 

(2.30) 

(2.31) 

(2.32) 


the  following  expression  for  the  bubble  scattering  cross  section  is  obtained: 
_ 7cQ'^p^(3(Y+1)(o§Q--a)f--o)|+coiO)2)^Pi^P2^ _ 

f^eqP?pi!(®on^^^r+5Qi^'^I(cooi^-(of)^+6fo)fJ(a)o2^-o)i)^+62a)^Pf-^+P2^ 

(2.33) 
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3.  DIFFERENCE  FREQUENCY  SCATTERING  CROSS  SECTION  FROM  A 
DISTRIBUTION  OF  BUBBLES  IN  AN  ACOUSTIC  FLUID 


Of  interest  is  the  scattering  cross  section  per  unit  volume  of  a  sediment 
that  contains  trapped  bubbles.  If  we  make  a  single  scattering  assumption,  based 
on  resonance  scattering  as  a  dominant  effect, 21  we  can  write  the  volume 
scattering  cross  section  in  the  following  form: 


P=|  on(RoJclRe, 


(3.1) 


where  Req  is  the  equilibrium  bubble  radius  and  n(Req)  is  the  bubble  size  density 
function,  defined  as 


ri(Req)  - 


dN(Req) 

dReq 


(3.2) 


where  N(Req)  is  the  number  of  bubbles  per  unit  volume  with  radii  less  than  Req. 

A  combination  of  Eqs.  (3.1)  and  (3.2)  results  in  an  integral  that  is  very 
difficult  to  solve  analytically.  A  technique  originally  employed  by  Wildt22  yields  an 
approximate  expression,  based  on  the  following  assumptions: 

1 .  Most  of  the  scattering  is  from  bubbles  with  radii  near  the  resonance  radius. 

2.  The  bubble  size  density  function  is  approximately  constant  in  each  range 
of  bubble  size  that  contributes. 

3.  The  damping  constant  d  is  approximately  constant  in  each  range  of  bubble 
size  that  contributes. 

Medwin's23  expression  for  the  resonance  frequency  fo  in  terms  of  the  equilibrium 
bubble  radius  is 
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(3.3) 


where  y  is  the  ratio  of  specific  heats,  Peq  is  the  ambient  pressure,  p  is  the 
ambient  density,  and  b  and  p  are  quantities  that  account  for  surface  tension  and 
thermal  conductivity. 

By  inversion,  the  resonance  radius  Rr  is  given  in  terms  of  applied  frequency  f: 


2itf 


By  combination  of  Eqs.  (3.3)  and  (3.4), 

h  _  Rq 

f  Req 


Therefore, 


*^01  _  i^01  0^02  _  i^02 

Req  ’  “2  Req  ■  ^  ~  Req 


(3.4) 


(3.5) 


(3.6) 


where  Roi  is  the  resonance  radius  of  a  bubble  driven  at  angular  frequency  coi. 
By  combining  Eqs.  (3.6)  and  (2.33), 


^pQ[3(Y+1)co§-o)f- 

-0)? 

J  — 

a)iO)2]^Pl^P2 

J 

[<  J 

2 

+  82 

fRoa 

ilReq  J 

2 

+  8^ 

M  ^2. 

(3.7) 

By  inspection,  it  is  clear  that  Eq.  (3.7)  will  have  maxima  when  R=Roi, 
R=Ro2>  R=Ron.  The  height  of  these  maxima  will  be  determined  by  the  values  of 
the  damping  constants  5i,  82,  and  Sq.  Assumption  (1)  is  a  statement  that  the 


12 


maxima  are  high  enough  that  most  of  the  contribution  to  the  integral  in  Eq.  (3.1) 
comes  from  bubbles  with  radii  very  close  to  Roi,  R02.  or  Ron.: 


^Roi+e 

P  =  Pi+P2+Pn  =  on(Req)dReq  +  an(Req)dReq  +  I  an(Req)dReq  •  (3.8) 

jRoi-e  ^  •/R02-«  •'Rocr« 

Consider  pn,  which  is  the  contribution  from  bubbles  with  radius  close  to  Ron.: 

^Roo+e 

Pn=  on(Req)dReq 

jRonre 


^R(Xjfe 


2d2 


P^PfB 


p?pio>t4p?5i+p||i 

l  ^2J 


[3(y+1)o)§  -  0)f  -  coi  -  a)i(oJn(Req)dReq 


r2 

Med 


fkifr® 


Rp 


1  +52 


^eq 


(3.9) 


According  to  assumption  (2),  the  bubble  size  density  function 
n(Req)=n(Ron)  is  constant  within  the  narrow  range  of  bubble  radii  about  Ron- 
According  to  assumption  (1)  and  Eq.  (3.6),  the  resonance  frequency  can  be 
approximated  by  coo=^i.  The  second  and  third  factors  in  the  denominator  of 
Eq.  (3.9)  don't  change  significantly  across  the  bubble  radius  range,  and  can  be 
approximated  as  constants  with  Req=Ron-  Hence, 


pQ= 


7tp2 

CM 

3 

T— 

3 

1 

CM  CM 

3 

1 

CMt- 

3 

1 

CM 

V?Pin(Ron) 

2  2  4  4 

PlP2“>2 

p2^n , p2^n 

Pl-^+P2^ 

(  R2  1 

M0I  1 

r2 

V^on  J 

2  1 

+8? 

f  p2  ^ 

"02 
r2 

V^on  J 

2 

+5i 

X 


^Rotl+i 


dR 


eq 


R 


eq 


HoQ-e 


R  2 
“eq 


-1 


+5 


(3.10) 
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To  solve  the  above  integral,  consider  the  variable  transformation 


q= 


^0£2_i 

Req 


dq  =  ' 


R 


on 


R 


dR 


eq 


eq 


(3.11) 


The  integral  in  Eq.  (3.10)  is  then  given  by 


^RoQ+e 


dR 


eq 


R 


eq 


RoQ-e 


R 


OQ 


•eq 


+52 


-q2 


-dq 


'on 


qi 


+8 


n 


(3.12) 


where 


qi 


Ron  1 
(Ro«-e)  ' 


q2  = 


Ron  .| 
(Ron+s) 


(3.13) 


Since,  according  to  assumption  (1),  the  contribution  to  the  integral  is  from  a 
narrow  range  of  bubble  radii  near  Req=Ron.  the  variable  q  is  always  small. 
Hence  the  following  approximation  can  be  made: 

(q(q+2))^  =  4q2  (3.14) 

Taking  advantage  of  this  approximation,  the  integral  is 


where  the  limits  of  integration  have  been  expanded  out  to  infinity.  This  is 
reasonable  since  the  contribution  to  the  integral  is  small  outside  the  original  limits 
qi  and  q2.  The  resulting  definite  integral  is  tabulated:24 
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I  -  1  71 


(3.16) 


Combine  Eqs.  (3.10),  (3.12),  and  (3.16)  to  get  an  expression  for 


Pi  and  P2  can  be  obtained  by  repeating  the  procedure  of  Eqs.  (3.9)  -  (3.17) 
above,  over  bubble  size  ranges  about  R01  and  R02.  The  results  are 


This  page  intentionally  left  blank. 
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4.  DIFFERENCE  FREQUENCY  SCATTERING  STRENGTH  OF 

A  GASSY  SEAFLOOR 

In  Ref.  6  an  expression  was  developed  for  the  linear  backscattering 
strength  of  a  sediment  interface,  due  to  trapped  bubbles  within  the  volume  below 
the  interface.  In  this  section,  the  same  treatment  will  be  followed  in  developing  a 
difference  frequency  scattering  strength. 

In  Section  4.1  an  expression  is  derived  for  the  scattered  pressure  from  an 
element  of  sediment  volume,  based  on  the  Biot  theory  and  the  principle  of 
acoustic  reciprocity.  In  Section  4.2,  an  expression  is  obtained  for  the  effective 
scattering  cross  section  per  unit  area  of  sediment  interface. 

4.1  SCATTERED  PRESSURE  FROM  AN  ELEMENT  OF  SEDIMENT 

VOLUME  BY  RECIPROCITY 

Consider  an  acoustic  source  in  the  water  column  and  an  element  of 
sediment  volume  that  scatters  sound,  as  in  Fig.  4.1.  The  source  and  scatterer 
are  both  assumed  to  be  much  smaller  than  the  acoustic  wavelength.  One  can 
surround  each  with  virtual  spheres  that  are  large  in  comparison  with  the 
wavelength,  as  in  Fig.  4.2.  The  spheres  surrounding  source  and  scatterer  will  be 
called  the  “source  sphere"  and  the  "scatterer  sphere",  respectively. 

The  source  sphere  starts  with  a  surface  velocity  Vq  ,  which  generates  a 
pressure  P  at  the  scatterer.  The  scatterer  sphere  responds  with  a  surface 
velocity  v,  which  in  turn  induces  a  scattered  pressure  Py  back  at  the  source. 
The  pressures  and  surface  velocities  are  related  by  the  principle  of  acoustic 
reciprocity  which  states  that,  in  a  linear  medium,  the  source  and  scatterer  can  be 
swapped  with  no  change  in  the  relationship  between  transmit  and  receive 
signals.  For  linear  scattering  problems  this  swapping  of  positions  can  be 
interpreted  to  represent  the  backscatter  case.  In  terms  of  pressures  and 
velocities,  reciprocity  states: 


H)"  V 


(4.1) 
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Source 


Water 

Sediment 


Scattering  element  dxdydz 


Figure  4.1 

Backscatter  from  an  element  of  sediment  volume. 
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The  surface  velocities  can  be  related  to  the  local  acoustic  pressures  and 
impedances: 


II 

(4.2) 

eP 

1#  ^  S£- 

^  Z  ' 

(4.3) 

where  Zq  and  Z  are  the  acoustic  impedances  in  the  fluid  above  and  in  the 
sediment,  e  is  a  transfer  function  from  incident  pressure  P  to  scattered  pressure 
eP  at  the  surface  of  the  scatterer  sphere.  It  includes  the  combined  effects  of 
scattering  at  the  element  dxdydz  and  propagation  from  the  scatterer  to  the 
surface  of  the  surrounding  virtual  sphere.  As  illustrated  in  Fig.  4.2,  this 
propagation  takes  place  as  if  it  were  in  the  water  column.  If  the  water  column's 
attenuation  is  neglected,  the  average  square  magnitude  of  e  is: 

=  ,  (4.4) 

s 


where  rg  is  the  radius  of  the  virtual  spheres  and  Obv  is  the  scattering  cross 
section,  defined  as  the  ratio  of  scattered  power  to  incident  intensity.  By 
combining  Eqs.  (4.1)  -  (4.4),  an  expression  is  obtained  for  the  square  magnitude 
of  the  pressure  returned  to  the  projector: 


P 

Pq/Zo 


2gbv 

47irs2 


jP/Zl^dxdydz 


(4.5) 


Since  the  sediment  under  consideration  is  a  Biot  medium  supporting  fast  and 
slow  waves,  the  acoustic  intensity  returned  to  the  projector  has  two  components: 


Pf+Ps 

H  %vf 

Pf 

%vs 

Ps 

J  Pq/Zo 

[Anri 

Zf 

Zg 

|2\ 

dxdydz 


(4.6) 


where  Obvf  and  Obvs  are  the  scattering  cross  sections,  per  unit  sediment  volume, 
associated  with  the  fast  and  slow  waves.  The  quantities  Zf,  Zg,  Pf,  and  Pg  are. 
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yo.Po 


(a)  Acoustic  propagation  of  incident  sound.  Virtual  sphere  surrounding 
the  source  has  surface  velocity  vq  and  pressure  po-  v-j  and  Pi  are  the 
pore  fluid  velocity  and  pressure  induced  at  the  scatterer. 


(b)  Acoustic  propagation  of  backscattered  sound.  Virtual  sphere  surrounding 
the  scatterer  has  surface  velocity  vi'  and  pressure  pi'.  P2  is  the 
backscattered  pressure  induced  at  the  source. 


Figure  4.2 

Calculation  of  backscatter  by  acoustic  reciprocity. 
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respectively,  the  fast  and  slow  wave  acoustic  impedances  and  the  acoustic 
pressures  carried  by  the  fast  and  slow  wave.  The  quantity  inside  the 
parentheses  is  the  square  magnitude  of  the  scatterer's  surface  velocity.  The 

factor  converts  this  to  the  backscattered  pressure  magnitude  squared. 


Pq/Zo 


Equation  (4.6)  involves  linear  scattering,  where  all  propagation  and 
scattering  occurs  at  a  single  defined  frequency.  In  the  case  of  parametric 
scattering  the  incident  sound  upon  the  scatterer  is  at  the  primary  frequencies, 
whereas  the  sound  propagates  away  from  the  scatterer  at  the  difference 
frequency.  Care  must  be  taken  in  specifying  the  frequencies  at  which  the 
pressures  and  impedances  in  Eq.  (4.6)  are  defined.  The  square  magnitude  of 
backscattered  pressure  at  the  difference  frequency  is 


dxdydz 


(4.7) 


where  pf  and  Ps  are  the  difference  frequency  scattering  cross  sections,  according 
to  Eq.  (4.1),  for  fast  and  slow  waves,  respectively.  Pit,  and  P2f  are  fast  wave 
acoustic  pressures  incident  upon  the  scattering  element  dxdydz  at  the  primary 
frequencies  coi  and  C02.  Pis  and  Pas  are  slow  wave  acoustic  pressures  incident 
upon  the  element.  Zif,  Zat,  Zis,  and  Zas  are  the  corresponding  partial  acoustic 
impedances. 

The  quantity  inside  the  brackets  is  the  Q  component  of  the  scattering 
sphere's  surface  velocity.  The  factor  ^  converts  this  to  the  backscattered 
pressure  returned  to  the  projector.  According  to  the  principle  of  acoustic 
reciprocity,  this  is  given  by 


Pfif+Pfis 


Pno^Zoo 


(4.8) 


where  Pm  and  Pqs  are  the  fast  and  slow  acoustic  pressures  that  would  be 
induced  at  frequency  Q  upon  the  scattering  element,  if  the  fluid  surrounding  the 
source  were  driven  at  frequency  Q  with  velocity  Pqo/  Zno- 
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4.2  DIFFERENCE  FREQUENCY  BACKSCATTERING  STRENGTH  OF 
SEDIMENT  INTERFACE 

In  the  following,  the  difference  frequency  backscattering  strength  of  the 
sediment  interface  is  defined.  The  definition  is  similar  to  that  for  linear  scattering 
from  bubbles.25  As  illustrated  in  Fig.  4.3,  the  effective  surface  backscattered 
pressure  Pg  per  unit  area  dxdy  of  sediment  interface  is  a  sum  of  contributions 
from  all  volume  scattering  elements  below  the  interface: 


<|p.R=j 


f 

IPtf'jPzf' 

'U  I 

Pl8 

2 

^23 

J. 

47rr|\ 

Zii  IZj, 

i  4jtr|l 

Zl8 

+ 

Zzs 

jj 

dz  .  (4.9) 


The  parametric  backscattering  strength  is  given  in  terms  of  the  pressure  Pjnc 
incident  upon  an  interface  element  dxdy  and  the  scattered  pressure  Psi  at  unit 
distance  rim  from  dxdy: 


,  (4.10) 

where  r  is  the  distance  between  the  source  in  the  water  column  and  the 
scattering  element  dxdy,  and  a  is  the  absorption  in  the  water  column.  The 
backscattering  strength  of  the  interface  is  defined  as 


BS=10log^^lS 

IPincP 


(4.11) 


Upon  substituting  Eqs.  (4.9)  and  (4.10)  into  Eq.  (4.1 1), 


BS  =  10log 


■r 

Pof+Pfls 

'f  pf  ( 

’’n  rj  Pa  n  4.  Ps  ( 

Pis 

2 

Pzs 

"11 

Pflc/Zoo 

47trs2i 

Zifn^lJ''47rr|l 

Zl8 

+ 

Z2S 

) 

fc)2 


mncP 


(4.12) 


22 


Projector 


Figure  4.3 

Effective  interface  backscatter  due  to  volume  scattering. 
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If  the  absorption  in  the  water  column  is  neglected, 


BS  =  10log 


Pof+Pns 

2 

Pf  / 

Pis 

2 

P2S 

Poi^^no 

4jtr|l 

J2„  IZa 

JZis 

*r 

Z2S 

|Pll+P21P 


,(4.13) 


where  Pii  and  Pr,  are  the  pressures  of  the  incident  sound  at  the  primary 
frequencies  ©i  and  ©2-  Poi  =  PnoCrim/r)  is  the  incident  acoustic  pressure  at  the 
difference  frequency  Q.  when  Pqo  is  generated  at  the  source  sphere. 

Figure  4.4(a)  is  a  plot  of  the  difference  frequency  backscattering  strength 
modeled  by  Eq.  (4.12)  for  a  typical  coarse  sandy  sediment.  The  sediment 
parameters  are  listed  in  Table  4.1.  The  backscattering  strength  in  this  figure 
increases  with  frequency.  This  is  consistent  with  expectations,  based  on  the 
bubble  size  density  function,  which  was  derived  from  the  measured  grain  size 
distribution  function.  Over  the  frequency  range  of  Fig.  4.4(a),  the  bubble  size 
density  function  increases  as  the  bubble  radius  decreases.  Since  the  model  is 
based  on  resonance  scattering  from  bubbles,  an  increase  in  the  scattering 
strength  with  frequency  is  expected. 

The  sediment  parameters  are  listed  in  Table  4.1.  The  backscattering 
strength  in  this  figure  increases  with  frequency  because  of  the  bubble  size 
density  function,  which  decreases  with  bubble  size.  Since  the  model  is  based  on 
resonance  scattering  from  bubbles,  this  means  that  the  backscattering  strength 
should  increase  with  frequency,  which  is  consistent  with  the  model's  behavior. 

In  Fig.  4.4(b),  the  dependence  of  the  backscattering  strength  on  the 
incident  primary  pressure  is  plotted.  The  two  primary  pressures  are  assumed  in 
this  figure  to  have  the  same  amplitude.  The  backscattering  strength  increases 
linearly  with  the  primary  signal  pressure. 
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(a)  Backscattering  strength  versus  difference  frequency  for  source  level  of 
200  dB  re  1  jxPa  and  center  frequency  of  42  kHz,  at  normal  incidence. 


(b)  Backscattering  strength  versus  source  level  (at  primary  frequencies)  for  a 
difference  frequency  of  10  kHz  and  center  frequency  of  42  kHz,  at  normal 
incidence. 


Figure  4.4 

Predicted  parametric  backscattering  strengths  over  a  gassy  sand 
with  geoacoustic  properties  given  in  Table  4.1. 
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Table  4.1 

Geoacoustic  input  parameters  for  parametric  scattering 
strengths  of  Fig.  4.4. 


Parameter 

Units 

Fluid  density 

1000 

Fluid  bulk  modulus 

2.25x1  o9 

Porosity 

0.4 

Grain  density 

(kg/m^) 

2650 

Mean  grain  diameter 

(<!))* 

1.0 

Standard  deviation 

(<!>)* 

1.0 

Pore  size  parameter 

(m) 

1.796x10-4 

Viscosity 

(kg/m-s) 

1.0x10-3 

Permeability 

(m2) 

6.45x1 0-^O 

Virtual  mass  parameter 

1.75 

Grain  bulk  modulus 

(Pa) 

7.0x1 09 

Frame  shear  modulus 

(Pa) 

2.61  x107 

Shear  log  decrement 

0.15 

Frame  bulk  modulus 

(Pa) 

5.3x1 09 

Bulk  log  decrement 

0.15 

Gas  bulk  modulus 

(Pa) 

2.48x1  o3 

Gas  density 

(kg/m3) 

1.22 

Gas  heat  conductivity 

(cal/m-s-®C) 

5.6x10-3 

Gas  spec,  heat  (const  press) 

(cal/kg) 

240 

Gas  specific  heat  ratio,  Cp/Cy 

1.4 

Bubble  surface  tension 

(N/m2) 

0.075 

Bubble/pore  volume  ratio 

0.625 

Gas  content 

1.0x10-3 

*  (j)  =  -log2(grain  diameter  in  millimeters) 


26 


5.  CONCLUSIONS 


A  model  for  parametric  backscatter  from  trapped  bubbles  in  sandy 
sediments  has  been  developed.  It  is  based  on  the  nonlinear  behavior  of  a  gas 
bubble  in  an  acoustic  fluid,  with  modifications  that  allow  the  treatment  of  bubbles 
in  a  poroelastic  medium.  These  modifications  consist  of  assigning  the  pore  fluid 
effective  densities  and  acoustic  impedances  that  characterize  the  Biot  fast  and 
slow  waves. 

The  model's  predictions  are  consistent  with  expectations.  The 
backscattering  strength  increases  linearly  with  source  strength,  which  is  identical 
to  the  case  for  parametric  scattering  from  bubbles  in  an  acoustic  fluid.  The 
increase  in  parametric  backscattering  strength  with  frequency  is  consistent  with 
the  bubble  size  distribution  inferred  from  the  input  grain  size  distribution. 

Similar  parametric  scattering  models  have  been  developed  for  bubbles  in 
water, 26  and  these  compare  reasonably  well  with  experimental  measurements. 
This  work  allows  an  application  of  this  modeling  technique  to  gas  bubbles  that 
might  be  trapped  within  sediment.  At  present,  however,  the  authors  are  aware  of 
no  experimental  measurements  of  parametric  scattering  strengths  from 
sediments  where  the  gas  bubble  size  distributions  are  known.  Such  data  are 
needed,  when  they  become  available,  to  test  and  further  develop  this  model. 
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APPENDIX  A 

DERIVATION  OF  NONLINEAR  EQUATION  OF  MOTION 
FOR  A  SINGLE  BUBBLE 
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The  following  derivation  follows  the  treatments  of  Noltinigk  and  Nepiras27 
and  Zabolotskaya  and  Soluyan .28  Their  papers  consider  the  volume  oscillations 
of  a  bubble  surrounded  by  water  and  insonified  with  a  parametric  source. 

Consider  a  bubble  surrounded  by  water.  We  are  interested  in  obtaining  a 
relationship  between  the  radial  velocity  of  the  bubble's  surface  and  the  applied 
pressure.  Let  the  bubble  have  a  radius  R,  which  differs  from  its  equilibrium  radius 
Ro. 


The  kinetic  energy  of  the  liquid  surrounding  the  bubble  is 


In  an  incompressible  medium,  the  radial  velocity  of  the  fluid  surrounding  the 
bubble  would  have  a  phase  independent  of  the  distance  r  from  the  bubble's 
center.  The  radial  velocity  dr/dt  would  then  be  given  by 


dt "  r2  dt 


(A.2) 


In  the  case  of  a  compressible  medium,  dr/dt  will  lag  in  phase  behind  dr/dt: 


dr_  r2  dR  p-if%)r 
dfT^  dT®  ^  ’ 


(A.3) 


where  X  is  the  wavelength.  Insert  Eq.  (A.3)  into  Eq.  (A.1): 
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When  the  bubble  radius  R  is  small  in  comparison  with  the  wavelength  X 
Eq.  (A. 4)  can  be  approximated: 


KE  =  2.p(R^ff£| 

=  aipR^^f  •  (A-5) 

This  kinetic  energy  can  be  equated  to  the  work  done  by  the  gas  inside  the 
bubble: 


Mfr-£ivp»» 

=  f  (P.-PHto^Wr 

jRo 


(A.6) 


where  Vq  and  Rq  are  the  bubble's  equilibrium  volume  and  radius.  Pa  is  the 
pressure  inside  the  bubble  and  P  is  the  ambient  pressure.  Differentiating  both 
sides  with  respect  to  R  gives 


atp(3R^^)\R5i(^)^)  =  (P,-P)(4nR^  ■  (A'!’) 


Upon  rearranging,  the  well  known  Rayleigh  equation  of  motion  is  obtained: 


dt2  ^2ldtj 


dt 


(A.8) 


where  R  is  the  bubble  radius,  p  is  the  density  of  the  surrounding  fluid.  Pa  is  the 
gas  pressure  inside  the  bubble,  and  P  is  the  ambient  pressure  outside  the 
bubble.  In  terms  of  the  bubble  volume  V  =  4/3  7tR^,  this  equation  can  be 
rewritten: 


=  Pg-P  , 
D  ®  * 


(A.9) 
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where 


a  = 


P0{4^)2/3 


V3 


(A.  10) 


Upon  making  the  substitution  aV-''/3=po/(47iR)  and  rearranging,  Eq.  (A.9) 
becomes 


o  1 4jtRp  _  47iRp 

''  6  V  Po  Po 


(A.11) 


If  we  assume  that  bubble  oscillations  are  adiabatic,  the  relationship  between 
bubble  volume  and  internal  pressure  is  given  by 


■ 


(A.  12) 


where  y  is  the  ratio  of  specific  heats.  For  small  oscillations,  the  pressure,  bubble 
radius,  and  volume  can  be  expressed  as  equilibrium  values  with  small 
perturbations: 


P  =  Peq+P' 
R  =  Req  +  P 

V  =  Veq  +  V‘ 


The  time  derivative  of  Veq  is  zero;  therefore. 


V  =  -^V‘ 
^  dt'' 


V  =  ^V 
dt^ 


(A.  13) 
(A.  14) 
(A.  15) 


(A.  16) 
(A.  17) 
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Equation  (A.  12),  for  small  perturbations,  is  therefore 


P  =  P 
*  a  *  eq 


Veq+V 


To  second  order  in  (V'A/eq),  this  is 


p  =p 


Substituting  Eq.  (A.  19)  into  Eq.  (A.11), 


V-l^-iirRp  _Y-^  + 

ft  \/  On  *^eq  /\/  ^ 


y'  .  Y(Y<-I)f  V 


6  V  Po  T  'V„  2  V, 


(A.  18) 


(A.  19) 


=  ^(Peq-P)  ■  (A.20) 


Noting  that  the  equilibrium  volume  is  given  by  Vgq  =  4/37cReq®, 

= fa(p^-p 

®  ^  PoReq  PoR4  SltRe^  Po  ) 

Upon  making  the  following  substitutions. 


(A.21) 


2  _ 

“^'PoRlq 

(A.22) 

B=  ' 

8ltR5„ 

(A.23) 

4ltReq 

Po 

(A.24) 

a  =  3(Y+1)co§p  , 

(A.25) 
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Eq.  (A.21)  can  be  written  in  the  following  form: 


(A.26) 


If  we  express  R  and  V  in  terms  of  their  equilibrium  values  and  perturbations, 
Eq.  (A.26)  becomes 


V  - ^  U  kv  -  aV^ 

®^eqlVea+V  I  \  ^eq  /V 


^Hcoy-aV'"  +  ePj  =  0.  (A.27) 


Vfia  Rfia"^^  . 

Expanding  the  quantities  - ^  and  —  to  first  order  inV  /Vgq, 

V-a+V  '^eq 

\  J  V  J 

»  -  ^(l-V'/V„)  +  (l  +V'/3Ve<|)|o)2v'  -  aV'%  epj  =  0 .  (A.28) 


Upon  rearranging,  this  expression  can  be  written  in  the  following  form: 


+  |a)§V‘  -  aV’^  +  eP  j  + 


(''"e-le 


(A.29) 


Equation  {A.26)  can  be  rearranged  as  follows: 


cogv’-aV'  +eP'  =  -^^  + W 


(A.30) 


Expanding  Req/R  to  first  order  in  VWeq, 


0)2v‘-aV‘^  +  eP'=  1-4x7^ 


-V  + 


(A.31) 
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By  applying  Eq.  (A.31)  as  a  substitution  in  the  last  term  in  Eq.  (A.29),  the 
following  expression  is  obtained 


i 


1  V  ,1V 


6VeqVeq'3Veq 


\/2^ 

''^6  V 


'-5I 


V 


V 


V, 


eq. 


=  0 


(A.32) 


By  throwing  out  terms  of  order  (V'A/eq)^,  Zabolotskaya's  equation  of  motion  for  a 
bubble  is  obtained; 


V--^|v^  +  2VV'j  +  o)2v'-aV'^  =  -eP‘  .  (A.33) 

This  expression  neglects  the  possible  loss  of  energy  during  bubble  oscillation 
due  to  damping.  In  order  to  include  damping,  a  term  proportional  to  the  time  rate 
of  change  of  bubble  volume  must  be  included,  resulting  in  Zabolotskaya's 
damped  equation  of  motion  for  the  bubble: 

V  +  co2V-aV2-p^2\/V+V2j+fV  =  -eP  (A.34) 

The  arbitrarily  defined  quantity  f  in  Eq.  (A.34)  is  related  to  the  bubble's  damping 

constant,  which  is  defined  from  the  linear  part  of  the  bubble's  equation  of 

motion.  If  the  nonlinear  terms  in  Eq.  (A.34)  are  neglected,  the  remaining 
equation  is  that  of  a  damped  harmonic  oscillator, 

V  +  fV  +  (ogV  =  -eP,  (A.35) 

with  a  damping  constant  given  by  , 

5  =  -Q  ,  (A.36) 

where  co  is  the  angular  frequency.  In  terms  of  5,  the  resulting  equation  of  motion 
(A.34)  is 


36 


(A.37) 


V  +  coJV  -  oV^  -  pf2'?V  +V^]  +  =  -eP 

In  Zabolotskaya  and  Soloyan's  treatment  of  bubbles  in  an  acoustic  fluid,  the 
coefficients  coo,  a,  p,  and  5  are  invariant  with  respect  to  co.  More  generally,  as  is 
the  case  for  trapped  bubbles  in  pores,  these  quantities  must  be  allowed  to  vary 

with  CO. 
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APPENDIX  B 

RELATIONSHIP  BETWEEN  ACOUSTIC  PRESSURE 
AND  VOLUME  OSCILLATION  OF  A  MICROBUBBLE 
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Equation  (2.23)  is  a  relationship  between  the  pressures  applied  to  a 
microbubble  and  the  amplitude  of  its  resulting  volume  variations.  We  are 
interested  in  the  acoustic  pressure  radiated  in  the  bubble's  farfield.  To  that  end 
we  need  a  relationship  between  a  bubble's  volume  oscillations  and  the  acoustic 
pressure  generated  in  the  farfield.  In  developing  such  an  expression,  the 
following  assumptions  will  be  made. 

(1)  The  bubble  is  small  in  relation  to  the  acoustic  wavelength. 

(2)  The  signal  amplitudes  are  small  enough  that  the  linear  wave 
equation  applies. 

The  linear  wave  equation  in  a  fluid  medium  is  given  by 

=  ^  ,  (B.1) 

where  c  is  the  phase  velocity  of  a  compressional  wave  in  the  fluid.  For  a  pressure 
field  that  is  spherically  symmetric  about  a  source  at  the  origin,  there  is  only  radial 

dependence.  The  operator  in  this  case  is  given  by 

=  •  P-2) 

Upon  substitution  of  Eq.  (B.2)  into  (B.1)  and  rearranging,  the  wave  equation  can  be 
written  in  the  form 


a^(rp)  ^  1  8^(rp) 
dr^  c^  3t^ 


General  solutions  for  the  quantity  rp  can  be  expressed  as 

rp  =  fi(ct--r)+f2(ct+r)  , 


(B.3) 


(B.4) 
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where  the  first  term  represents  a  general  incoming  wave  and  the  second  term  a 
general  outgoing  wave.  If  we  consider  outgoing  harmonic  waves  only,  the  pressure 
field  can  be  expressed  as 


(B.5) 


where  A  is  an  arbitrary  complex  amplitude,  k=co/c+ia  is  the  acoustic  wavenumber,  a 
is  the  absorption  in  Np/m,  and  (o=27cf  is  the  angular  frequency.  For  small  amplitude 

signals  the  fluid  velocity  is  related  to  the  acoustic  pressure  according  to  Euler's 
equation 


Po^  =  -Vp  .  (B.6) 

where  po  is  the  average  density  of  the  medium.  Upon  combination  of  Eqs.  (A.5) 
and  (A.6),  a  relation  between  acoustic  pressure  and  fluid  velocity  is  obtained: 


I 

krjPoC  ■ 


(B.7) 


According  to  assumption  (1),  the  quantity  ka  is  small,  where  a  is  the  bubble 
radius.  Therefore  the  second  term  in  Eq.  (B.7)  dominates  at  r=a: 

P(a) 

kapoC 

The  acoustic  pressure  at  r  is  inversely  proportional  to  r: 

|pW|  =  f|p{a)|  ■ 

By  combining  Eqs.  (B.8)  and  (B.9)  and  substituting  lkl=co/c, 

I  I  a)a2po|v(a)| 

iPWh - r-^  • 


(B.9) 


(B.10) 
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The  rate  of  change  of  the  bubble's  volume  perturbation  V  is  related  to  its  surface 
velocity  v(a): 


^  =  47ta2v(a).  (B.11) 

If  the  oscillations  are  harmonic,  the  bubble  volume  rate  dV/dt  is  reiated  to  its 
volume  perturbation  V: 


dt 


=  a)V  , 


(B.12) 


where  co  is  the  angular  frequency.  Upon  combining  Eqs.  (B.10),  (B.1 1),  and  (B.12), 


Q>^Po|V| 

47tr 


(B.1 3) 


This  is  Eq.  (18)  in  Zabolotskaya  and  Sutin's  paper,^^  originally  derived  by  Landau 
and  Lifshits.34 
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